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Abstract 

An analysis is given of particlelike nonlinear bound states in the Newtonian limit of the coupled 
Einstein-Dirac system introduced by Finster, SmoUer and Yau. A proof is given of existence of these 
bound states in the almost Newtonianian regime, and it is proved that they may be approximated by 
the energy minimizing solution of the Newton-Schrodinger system obtained by Lieb. 

1 Introduction 

Since the gravitational interaction is attractive, there is the possibility that theories involving fields coupled 
to gravity will possess bound states held together by this force. At the level of Newtonian gravity an example 
is provided by the Newton-Schroodinger system, in which there are bound states which can be obtained by 
minimizing a nonlocal energy functional (|10|): see also more recent discussions in |H1 1121 [5]. In general 
relativity, it was first established that this phenomenon occurs for the Yang-Mills-Einstein equations (see 
[U [13] for the original work, and [3] for some more recent developments and further references) and then 
for various other systems involving Dirac fields ([Sl[7])- In this article we discuss the bound states in the 
Einstein-Dirac system which were studied numerically in [5] , and prove that in the almost Newtonian regime 
they may be rigorously approximated by the bound states in the Newton-Schrodinger system, and as a 
consequence give an existence proof. 

The starting point is the action functional 

S=^jRdl^9 + jZ^^j^A{p-'^)^Adn, (1.1) 

describing the interaction of two Dirac spinor fields ^'i and \l/2 with a gravitational metric g, whose scalar 
curvature is R and whose volume element is d^g. We use J?) for the Dirac operator derived from g, with 
associated 7 matrices, as explained in ^ Section II]. The Euler-Lagrange equations are 

„ 1 „ 87rG„ ,^ 7nc^ 

Rab - ^Rgab = ^Tab, {P - ^)*A = 1.2) 

where Rab is the Ricci curvature and Tab = \^ A{'i^adb + i^bda)'^ a\ is the energy-momentum tensor; 

space-time indices a,b take values in {0,1,2,3}. The system (|1.2p is precisely that studied in 0, except 
that the dimensional constants G, h, c have been reinstated. It is always to be understood that /i, m, G are 
fixed positive numbers, while the speed of light c is a positive number taking large values, i.e. we will study 
the Newtonian (or non-relativistic) limit in which c — > -l-oo. The reason there are two Dirac spinor fields 
in the model is that this allows for spherical symmetry: indeed, in the article 6 a spherically symmetric 
ansatz was introduced, and the corresponding system of ordinary differential equations (ODEs) for static 
spherically symmetric solutions was derived. Numerical evidence was presented for the existence of solutions 
to this system of ODEs, which give particlelike, or solitonic, solutions of p.2p corresponding to nonlinear 
bound states of two fermions; also the linear stability of these solutions was investigated. In this note we 



will show that in the Newtonian limit c +00 these solutions of the system (11.21) can be approximated by 
solutions of the Newton-Schrodinger system 

ih^ ^--A^P + u^:, -Au^-kIV-P (1.3) 
ot 2m 

for K = 87rGm. In fact the rigorous analysis of the limit we provide yields a proof of existence for the bound 
state solutions of (|1.2p introduced in [B], which does not seem to have appeared previously in the literature. 
The type of solution of Ijl.Sp is one in which the time dependence is uniform phase rotation at frequency j- : 
explicitly -ip — e~''^*(^(x), and ip solves 

^2 

Tjf — —- — A(p + mucp, — Au ^ —8TrGm\(p\ . (1-4) 
2m 

The existence of such nonlinear bound state solutions to the Newton-Schrodinger system was proved in |10j 
by variational methods. Pseudo-relativistic generalizations of Lieb's solutions have been given in [9], and we 
will make use of a result from this article on the non-degeneracy of the linearization of (|1.4p : see lemma S] 
and the subsequent remark. 

1.1 Spherical symmetry 

Using the ansatz for bound states with spherical symmetry introduced in [5], we search for solutions with 
metric 

g = c^e^^dt^ - e^^dr^ - r'^dif (1.5) 

and spinor fields of the form 

^ = e'-Z^e— ' ( V *2 = e^/^e— * ( .^^^^ \ (1.6) 



where 















62 = ( 


n 



1 ^ 

and fj' = — > 
J* z ^ 



where cr* are the Pauli matrices. In (|1.6p v, A, $1, $2 all depend on r only, and 

uh — mc^ + rj 

with -q as in (|1.4p . 

Remark 1 The angular dependence in (|1.6p is taken over directly from that displayed by the ground state 
Dirac wave functions for the relativistic hydrogen atom. Recall from USH 1161/ that in problems involving the 
Dirac equation with spherical symmetry there are three commuting operators: J^^Js and K (respectively 
the total angular momentum squared, the third component of the total angular momentum, and the spin- 
orbit coupling operator defined in UtA Equation 3.275]) which also commute with the Hamiltonian. The 
simultaneous eigenspaces of these three operators (corresponding to eigenvalues conventionally written j{j + 
l)fi? , j^h, nh) are two dimensional and give the decomposition into irreducible subspaces. The possible values 
of j are in the sei {i, |, . . . }, those of js are in the set {—j, . . . , +j} andK = ±{j + ^). The ground state 
corresponds to the choice j ~ i, and the wave functions V^i (resp. '^2) correspond to js = ^ (resp. J3 = —■^) 
and K — —1. The value of k is not important, but what is crucial is that the states ^i,'^2 have opposite 
values of j^, which ensures that the energy momentum tensor Tab is consistent with the spherically symmetric 
metric in (|1.5p ." see J^, Setion IV]. 
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Substitution of (|1.5p - (|1.6p into (|1.2p leads, as in [SI Section IV], to the following system of equations: 



^ ^j*!"*^ 

c n or r 



)$i-e-^— ^-^ ^*2-0, (1.7) 



)^2 + e-^—^ + ^ ^$1=0 (1.8) 

' or 



to be solved coupled to 



e 



-2A 



(2rA' - 1) + 1 = SvrG^p, (1.9) 
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e-2^(2ri^' + l)-l = 87rG-rp, (1.10) 

where p = p(i^, $i, $2, ??, c) and p = p{v, $1, $2, ''7, c) are given by 

p = 2(m + c-2?7)e"2''($? + $^) (1.11) 
]5 = -2(m + c"\)e"2.^($? + $^) (1.12) 

+ 4— e""*!*, + 2mc2e"''($?-$^) (1.13) 
r 

= -277e-^''($? + $^) - 2mc2e-''$2(g-i. _ ^^4) 

he 

- 2mc2e"''$^ +4— $i$2- (1-15) 
r 

1.2 Newtonian limit 

We will solve the system (|1.7l) - p.l0p in the almost Newtonian regime. To be precise, we treat e = as 
a small positive parameter, and show that the system (|1.7p - (|1.10p can be regarded as a perturbation of the 
Newton-Schrodinger system 

The spherically symmetric Einstein- Dirac system (ll.7p - (|1.10p admits nonlinear bound state solu- 
tions (A', i^', $1, <i>|) for small positive e, which can he approximated (in a strong weighted norm) 
by the bound state solution ip of the Newton-Schrodinger system which minimizes the energy 
(|2.25p .' in particular, (<i>|, $|) converges uniformly to {ip, 0) as e ^ 0. 

To formulate more precisely and prove this we carry out a rescaling of the dependent variables, and also 
make an adjustment to facilitate the handling of the ADM mass. The precise statement appears in theorem 
[9l and asserts the validity of the Newtonian approximation in norms stronger than with exponential 
weights in the (massive) Dirac fields $i,$2, and polynomial weights for the (massless) metric components 
X,v: see p.43p . The bound state of the Newton-Schrodinger system determines space and time scales, and 
the requirement that e be small (or, equivalently, that the speed of light c be large) should be regarded as 
being relative to the scale so determined. 

1.3 Auxiliary conditions 

We will work with boundary conditions corresponding to an asymptotically Minkowskian metric, i.e. 

lim e^^^"^) = 1 = lim e^''^''), 

r— *+oo r — >-+oo 



and also require finite ADM mass, so that 



hm ^(e^^W -1) = £2; 

r — *+oo 2 
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exists and is finite (for each positive e = c~^; the scaling factor is introduced for later convenience in 
analyzing the limit e — > 0). 

As r —> the solutions will satisfy e^^ — 1 + O(r^) and <i>2 = 0{r). The normalization condition 
II^A||L2(|j^j,o„^tQ„ti,') = 1, which ensures that the total probability equals one, can always be imposed on 
solutions by rescaling, see Section VI]. We will therefore not impose it throughout this paper, since a 
simple rescaling ensures that the solutions we obtain in theorem [S] do satisfy it. 

1.4 Rescaled variables 

We now introduce new variables Q,N,ip2 (in place, respectively, of A, i^, $2) which take into account both 
the auxiliary conditions and the expected behaviour in the Newtonian limit: 

e-2(i_e-2A)^2//o(r) + Q e"',^ = N, (1.16) 

with fo{r) = r^/(l + r)^, and 

^2 = eV'2. (1.17) 



Remark 2 The variable I determining the ADM mass is determined dynamically, as will become clear in 
the implicit function theorem set-up following p.36p . It turns out to be convenient for the analysis, however, 
to separate it off by introducing the function /o and requiring Q = 0(r"^) as r +00. Other choices of fo 
with the same asymptotic behaviour would be possible: the choice made here is just a simple function having 
the right asymptotics as r —> 00, and vanishing to 0{r^) as r — s- as is natural, see remark^ Notice that I 
is only uniquely determined by the requirement Q = 0{r^^) as r — > +00, and this condition is built into the 
function spaces used in the implicit function theorem. 

1.5 Formal consideration of the Newtonian limit 

Equations (fLg|) - prTU| ) become 

{rQY + 2lirfoy -SnGr^Pe^O, (1.18) 

l~e^2lfo + Q) 1- 62(2^/0 + Q) ^'-'"^ 

with p,{N, $1, V2, 7?) = pie^N, $1, eV'2, and p,{N, $1, V2, ??) = pie^N, $1, eV'2, while ^-i^ 

become 

(|-^)$.-^-^^2-F,^0, (1.20) 
h ' or r 

^^2 + ^-F2=0, (1.21) 
n or 



where 



and 



i^i - 2 ^$1 1.22 

+ ^(l-e-^^-)<i>r-e^(Q + 2^/o)(^ + i^2) 



F2 = e\Q + 2//o)(^ + -$1) - e2|e-^^V'2 + ^(1 - e-^^)^2. (1-23) 
^ or r ' n n 
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Note that in the Hmit e — > 0, the inhomogeneous terms Fi and F2 are formally O(e^), while 
= 2m|$i|2 + 0(e2) 



(1.24) 



r 



o(i), 



so that (|1.18p - (|1.19p imply that N becomes the Newtonian potential (i.e. the function u in (|1.4p ). and 
(eliminating -02 between (|1.20p - p.2ip ) the function $1 approaches (p, a solution of the Newton-Schrodinger 
system ()1.4p . So far this is completely formal, but it will be made into a rigorous statement in theorem [9] in 
section [3] In order to achieve this it is necessary to recall some facts about (II. 4p . 

2 The Newton-Schrodinger system 



We briefly summarize the approach to (|1.4p adopted in the article [10] so as to provide sufhcient details for 
our purposes. 

Theorem 3 ( |10j ) The nonlocal energy associated to (|1.4p 

^ I |V^(.)pdx -m^cf ( M^feM^dxd, (2.25) 

J J J \^ y\ 

admits a finite lower bound subject to the constraint of having J \(p{x)\'^dx — 1 fixed, and this lower bound 
is attained on a function which is unique up to translation. Further this minimizer is positive, spherically 
symmetric and a monotone non-increasing function of the radial coordinate satisfying \'p{r)\ < cie^'^^'' for 
some positive numbers c\,C2- 

We summarize the points of the proof from [lOj which will be needed below. 

1. The existence of a spherically symmetric minimizer of the nonlocal energy (|2.25p is proved by means of 
the Riesz rearrangement inequality, and a strict version of this inequality implies that any minimizer 
is spherically symmetric. The corresponding Euler-Lagrange equation is 

A^(a;) - 2m^G [ ^^^^dyifiix) = rjipix) (2.26) 



2m J \x — y\ 

where 77 < is the Lagrange multiplier. 

2. The relation between (|1.4p and (|2.26p follows quickly using the condition limi^jj^+oo |u(a;)| — and the 
formula for the solution of Poisson's equation —Au = f on R'^, namely: 

(-A-V)(x) = / TT^l^y- (2.27) 
J 4:TT\x-y\ 

3. For the case when f{x) — Kp{r), where k > is a constant and p is a function of the radial coordinate 
r = only, a result of Newton (jTT]) implies that (|2.27p can be rewritten: 

■ds 



(-A ^Kp){r)^-K.j (^^ -^'jp{s)s^ds + K J p{s)i 

= / ( )p{s)s'^ds + u{0). 

In r/ 
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Define K{r, s) = 87rs^(i- — i); this kernel is non-negative for < s < r. It is shown in 'TU] that any 
energy minimizing solution to (|1.4[) is radially symmetric, and so solves the equation 

E^ = -^{^ + ~y + ^'G Kir, sms)\^ds) ^ (2.28) 

where E — [rj — mu(0)). 

4. It is proved in [10] that all positive solutions of (|2.28|) can be obtained by a scaling of the unique 
positive solution of 

- + (i' K{r,sms)\^dsy = 0. (2.29) 

(The proof of uniqueness of positive solutions of (|2.29p is a crucial part of the article jTOj). Conversely, 
as long as £' > 0, there is a positive solution of (|2.28p . which can be obtained by scaling the unique 
positive solution of (|2.29p . 

Lemma 4 ([9]) Let L be the linear operator obtained by linearizing (|2.26|) ; 

L,(.)= (-I-A -,),(.) -2™^G / |fi^.,x(.)-4„.^G / i^j^.,,(.). 

Then L is a linear homeomorphism from H^^^ to L^^^ (where these are the radial parts of the usual Sobolev 
and Lebesgue spaces). 

Proof Recall that 77 < 0. It is proved in |9j that the operator L is self-adjoint on L^^^ with domain H'^^j^ 
and trivial kernel and a single negative eigenvalue with one dimensional eigenspace. (Going out of the radial 
sector, L has a three dimensional kernel associated with translation invariance, but this does not concern us 
in this article.) From this it follows that there exists c > such that H-Lxlli^ > c||x||l2, and that L maps 
^rad continuously onto i^ad- Furthermore the standard elliptic estimate gives ||L^^/||/f2 < C||/||l2 so that 
L is a linear homeomorphism H"^^^ —> L'^ad- ^ 

Remark 5 Curiously in the present general relativistic problem it is possible to use diffeomorphism invari- 
ance to circumvent the need for this lemma. Briefly, by rescaling time, it is possible to linearize the system 
with the value of the gravitational potential fixed at the origin r = rather than at infinity; for this lineariza- 
tion the corresponding nondegeneracy result is rather easy to prove. This was the approach taken in the first 
version of this article, and it has the advantage of not needing the more careful analysis of L given in f^. On 
the other hand it has the disadvantage of giving solutions which are not asymptotically Minkowskian ( without 
further co-ordinate changes), and this makes the physical picture less clear in various ways. Therefore it 
seems preferable to make use of the results of and obtain solutions which are immediately asymptotically 
Minkowskian. 

3 Theorem on the Newtonian hmit 
3.1 The analytic set-up 

We now introduce the spaces of functions with which we will be working. Let BC^ = BC^{[0,oo)) be the 
space of bounded continuous functions on the half line < r < 00 with H/Usc'' = sup(l + rY\f{r)\ < 00; 
notice that elements of BC^ satisfy / — 0{r^^) as r — > 00. Let i?C|, be the subspace of BC^ consisting of 
functions / vanishing to order S' > at the origin, in the sense that 

ll/llsc^., = supr-^' 1/(01+ sup(l + r)*'|/(r)| < 00. 

Further, let BC^'^ — BC^'^{[0, 00)) be the subspace of BC^ consisting of functions / which are also contin- 
uously differentiable with /' e BC^^^ , with the norm 

ll/llsci,. = sup(l + r)VMI + sup(l + r)i+*|/'(r)| < 00, 



6 



and, for S' > 1 let BC^f C BC^-^ be the subset of those vanishing at the origin to order 5', with norm 
ll/llsc^'"' ~ ll/llsci '5 + II/IIbc'', + II/'IIbc''+^ ^ ^rad ^'^d L^^^ be the radial parts of the usual 

Sobolev and Lebesgue spaces. We will use without comment the following inequalities: there exist positive 
constants Ci, C2, C3, C4, i?, such that 

sup |/(x)| <Ci||/||^2(R3) (3.30) 

|1/|1l6(R3) < C2||/||//1(r3) (3.31) 
ll//k|||L^(E3) <C3||/||h1(r3) (3.32) 

sup |a;/(x)| < C4||/||^i(i{3') for radia/functions / (3.33) 

aneR=':|x|>fl,. 

(see [TS], or [3]), and also the fact that H'^ad^'^) functions are continuous and 0(r~^) as r — > +00. The 
inequality p.32p is Hardy's inequality. The rad suffix will be ommitted when writing a norm. 

We also need the exponentially weighted Sobolev space H^'^'^^ defined, for s €E {0, 1,2,...} and (5 > as 
the subspace of L^(R'^) with the norm 



\f\\Hi^.^^= ||e*l-l9"/ll 



■m:|m[— 



finite; we write H^l^P for the radial part of H^^'^^ . Since / € Hj:^'^ immediately implies e'''^'/ € H^ady 
with a bound for the norm, we have by (|3.33p : 

sup |a;e*'''^'/(a;)| < C5||/||j:^{i,5} for radiaZ functions /. (3.34) 

xGR3:\x\>R, 

Introduce the Banach spaces 

X^Rx BCl^"" X BC^^^ X {hIIP n BC^'°) x (i/^^^^ n BCi''") (3.35) 
Y = BCl X BC' X {Hj:lf n BC) x {Hj:"J^ n BC) , (3.36) 

where for an intersection of two normed spaces we use the norm |j/||AnB = ll/IU + I!/I!-B- Now notice 
that (|1.18p - (|1.2ip define a system of equations .F(Z, Q, TV, $1, -02; e) = 0, where T : X x R ^ Y maps 
{I, Q, N, $1, tp2) G X and e G M (assumed small) to 

I (rQ) + 2/(r/o) - 8^Gr p„ 2rN - ~ —— 



l-e^2lfo + Q) l-e2(2Z/o + Q)' 
2m 94>i .77 aV'2 2 ' 

-^1p2 + ^ F2, (- '1P2-F1 

n or n n or r 



in Y. 



Remark 6 As promised in remark\^ the field Q is 0{r ^) at infinity (since we are solving for Q e BC'2'^ C 
BC^ ), so that the parameter I encoding the ADM mass is well defined, by integration of the first equation: 

2/= / SirGr^p.dr. 
Jo 

Remark 7 Notice that Q e BC^ , so that Q = 0{r'^) as r 0. To see that this is natural, integrate the 
first equation up to r, to deduce rQ + 2lrfo{r) — SttGs^ p^ds — 0{r^), and recall that we chose fo to be 
0(r2) asr-^O. 



7 



Remark 8 The fact that !F{ ■ ; e) maps X into Y (for small e > 0) can be read off from its definition and 
the formulae (|1.22[) - (jl.24p . by means of the inequalities (l3.30[) - p.34p above. For example consider the third 
component of J-': -02 cind are clearly in Hj:^'^^ H BC, so it remains to consider ¥2- Referring to ()1.23p . 
it follows from the fact that both Q and fa{r) are O(r^), that the first term is continuous at r = 0, while 
the same is obviously true for the second and third terms, since ■ip2 € BCl'° and N € BC^^^ . Furthermore, 
functions in H^^'^^ are continuous away from the origin and exponentially decreasing by (I3.34p . so that 
altogether F2 £ BC. Similarly F2 € H^^'P : the only term for which this is not immediately evident is 
but this is premultiplied by Q + 21 fo = O(r^) as r —> 0. These observations show that T is continuous from 
X X K ^ y (locally, near e = j; this will be strengthened to locally C^ in lemma 

3.2 The Newtonian limit and statement of the main theorem 

When e = the formal Newtonian hmit solution mentioned previously corresponds to the observation that 

T{2mG, 2rv! - 4mGfo, u, ip, O) = (3.37) 

2m 

where (p and u solve (|1.4p . To see this, observe that by (|1.24p and the preceding remarks, the equation 
J^(Z,(3,iV, $i,'02;O) = amounts to 

(rQ)' + 2/(r/o)' - 16rmTGr^\^if = 0, (3.38) 
2rN' ~ {21 fo + Q) = 0, (3.39) 
,77 mN ^ dib2 2 

U-^)^i--^--^2^0, (3.40) 
a n or r 

The first two equations imply — AA^ — — BtttoGI^iP, and substituting this into p.40p and then substituting 
for ip2 from (|3.4ip . implies that (<i>i,iV) solves (ll.4p . We choose, for e = 0, ($i,Af)to be the energy 
minimizing solution (</?, u) of (jl.4p described in section O The multipole expansion for N = u and the 
requirement that Q S BC^'^ imply that I = 2mG, since J \(p{x)\'^dx — 1. Collecting this together we write 
Stv — (2mG, 2ru' — AmGfo, u, ip, —^'p') and call this the Newtonian limit point in X. We can now state 
the main theorem: 

Theorem 9 There exists an interval (— ei, +ei) on which is defined a C^ curve e ^ 'Ef- = (Z*^, Q*^, A'''^, ip2) & 
X of solutions to the Einstein- Dirac system (jl.l8p - (|1.2ip such that ||S"^ — S^rUx = 0(e). More explicitly, 
l'^ — > 2mG and 



\Q' + AmCfo - 2ru'\\^ci.2 + \\N' - u\\bc^.i + ||$i - vWniJfnBC^." (3-42) 

h 

2^1- 



In terms of the original variables of the problem, we define a metric 

/^e-2e2-^dt2_g2A^^^2_^2^^2 

where u'^ — t^N'^ and (1 — e^^'^') = e^{2l'^fo{r) + Q'). Define also $| — e^'h then for e — c^^ small we have 
a solution (A',i^%$f,$|) to (frTjl - lLTOl) . and 

e-2||(l _ e-^^') - 2e^ru'\\^ci.2 + t-^\\v' ~ e^u\\Bc^.i (3.43) 
+ 11*1 - ip\\Hil-pnBC^^. + (^'^W2 + ^^v'\\Hll:PnBcl-° = 0[e). 
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Remark 10 We briefly discuss regularity of the solutions just obtained, omitting the e index for clarity. 
The definition of the last two factors in X ensures that $1, ■02 are for r > 0, and ip2 = 0{r) near r = 0; 
notice that (|1.2ip implies that = 0{r) as r — > also. It follows that the corresponding Dirac spinor 
given by (II. 6p is also on R"^. To see this, first oberve that the partial derivatives at r ^ of the function 
x'^il}2/r exist and given by 9j(x*V'2)(0) = (5ij'02(O)- -^"^ '^^■"^ calculate 

and notice that this is continuous everywhere, including at r ~ since 

ip'^iO) ^ lim lim 7/>^(r). 

r—>Q T r— (0 

Also the functions e^^ and N are functions for r > Q, whose derivatives have limit zero as r 0, and 
hence define functions on R"^ . Thus overall we have a classical solution of the equations of motion. 



3.3 Linearization and proof of main theorem 

Our aim is to use the implicit function theorem to obtain solutions of = for small e, thus proving theorem 
m To achieve this it is sufficient to check that T is locally and its partial derivative in the X direction, 
DiT , evaluated at the Newtonian limit is a bounded linear bijection (and hence a linear homcomorphism 
by the open mapping theorem). Let be the ball of unit radius centered at S in the Banach space X. 

Lemma 11 For sufficiently small the map T is in the neighbourhood i3i(Sjv) x (— e,,e») of the 
Newtonian limit point (SAr;0). The partial derivative DiJ^(Sjv;0) is the linear map X ^ Y which takes 
5 = {Sl,q,n,xi,X2) e X to 

[rq)' + 25l{rfo)' - 327rmGr2((^,xi), 2m' - 25//o - q, 

2m dxi /V mu^ 8x2 2 mn \ 

in Y . 

Proof This can mostly be read off by inspection. Consider for example the first component: {rQ)' +2l{r fo)' — 
SirGr^Pf. This will define a smooth map from X x K to BC^ if each of the three terms does so separately. 
Clearly Q [rQ)' is a continuous linear map BC\'^ BC^ , while the function {rfo{r))' lies in BC^ , so 
that I ^ l{rfo)' is continuous and linear from R to BC^ . Next, from the formula (|1.24p it is clear (using 
the inequalities (|3.30p - (|3. 341) 1 that € SC^, and that the corresponding map is (in fact smooth) from 
X X R to BC^ . The second component of T is handled similarly, after choosing e, to be sufficiently small 
that |1 — €^{2lfo + Q)! > ^ on Bi{En) x (— e*, e*). The third and fourth components can also be handled 
in a straightforward way, noting that the coefficients in the expressions ()1.22p - ()1.23p (i.e. the quantities 
multiplying ^i,ip2 and their derivatives) are smooth hmctions on BiCE.^) x (— e*,e*) taking their values in 
BC^'^; see also remark |51 □ 

Lemma 12 Let ip be an energy minimizing solution of p.4p as described in theorem\^ and let S^r be the 
corresponding Newtonian limit point satisfying !F{^m]Q) = as just described. Then the partial derivative 
in the X direction, Di!F{^m] 0), is a linear homcomorphism X ^ Y . 
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Proof To establish this we need to prove the unique solvabihty (with bounds for E X) of the equation 
-DiJF(Sjv; 0)(f) = y, for y = {a,(3,ji,j2) G Y. Thus we consider the system 

(rq)' + 25l{rf^y - 327rmGr2 xi) = a, (3.44) 
2rn - 251 -q^ (3.45) 

— X2 + ^=Ji, (3.46) 
,ri mu. dx2 2 mn , 

To start with, multiply the second equation by r, differentiate, add the first equation and then divide by 
2r^, leading to 

Now use the representation (|2.27p for n, and substitute the resulting formula for n into the equation obtained 
by eliminating X2 between the third and fourth equations: 



An+167rmG(v3,Xi> + — ttV^ = 0. (3.48) 



^2 ^ ^2 / 9 2\ 

Axi - VXi +'muxi + mrvf = ~hj2 - t~ + ~ -^i- 
Zm Zm \or r I 

This yields 

Lx\ = J 

where L is the operator in lemma |4] and J is the function 



J(r) = -h32 - 1^ + 

2m \ or r 



VTLLp 



(-A)-i( 



a 



(r/3)'. 



2r2 



In order to apply lemma [4] we need to check that J G i^ad- ^^i^ ^^^^ two terms clearly are, so we just 
need to show the same is true of the third term. Let tl).^^ be the homogeneous Sobolev space with norm 
11/11^1 = / l^/Pf^^^;, and let be the dual space. Recall that if^ad ^rad ^i^^ ^ continuous embedding, 
by Sobolev's inequality. Now ajr e so r~'^a £ H~^^ since, by Hardy's inequality. 



I / f^r^dr\ < \\a/rh4f/rU. < qi«|lBC^||/||^i. 



(Here the norm is understood to be with respect to the measure Anr^dr since all functions are radial.) 
Also r^'^{rpy e H^^^^ since 

I / f^^^^drl = I J r(3f'dr\ < ||/3/r|U. ||/'IIl^ < C||/3||bc^ 

Given this it follows from the Riesz representation theorem that (— A)^^ (^^^^?^^) G H^^^ C i^ad well- 
defined, and the bound 



holds. Therefore ip{—A) ^[ "'^r'^' ) G L^^^ (since iy9 is smooth and decays exponentially and so is in L^^^, 
and 5 = 5 + ^)- As a consequence there exists a unique xi G ff^ad such that Lxi = J satisfying the bound 

Ilxi||ff2 < c(^\\{a,(3)\\Bc^xBC^ + Hhh^ + lliillffi)- 
To derive the exponentially weighted bound, we apply lemma [T51 below to the equation: 

2m [ It? 2 

(-A + M + Vir))xi =.f = -f^ [-fij2 - + -n) - mmp 
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where M — —2rari/fi^ > and V = 2m^u/h'^. This imphes 



||Xl|l//{2.^} < c(^||Xi||h2 + ||/||if{0.^}) 



< c 



BC^xBC^ + \\:h\\Hi«-^i + lliillffd-''} + \\'n\\L'^\\e^''v\\L3 



for S small {S < \fM and less than the decay rate of in theorem [3]) . The third equation above then 



determines xi S -ff^ad'^^ uniquely: X2 = ^{ji - with the same bound for ||x2|iff{i,«} as for ||xi||//{2,5 



9X1 



We show that also X2 € -BC|'° below. First we must discuss q and n. Since we are attempting to solve 
for q = 0{r~^) as r ^ +oo, integration of the first equation gives 

/"OO 

2Sl= {32TTmGr^{Lp,xi) + a)dr 
Jo 



and 



qir) = 2Sl{- ~ /o(r)) - - I {327TmGs^cp{s), xi{s)) + a{s))ds 



-2dl(s/o(s))'+ (327rmGs2(^(s),Xi(s)) +a(s)) 
-2Sl{sfo{s)y + (327rmGs2(^(s), xi(s)) + a(s)) 



ds 



ds 



Since xi is bounded and continuous, the first line makes it clear that q — 0{r^^) as r ^ +oo, while the 
third implies that lim^^o exists (since a G i?G|), and hence that q G BC\'^ . Recalling that ip 
is exponentially decaying, and putting together the obvious estimates gives the bounds \5l\ < c(||a||BC'2 + 
llxillffs) and 119113(71 < c{\Sl\ + ||q;||bc| + llxillff^)- The similar bound for H^'Hscf then follows from (13.44^ . 
yielding the bound for ||(?||i,2;2 the norm on BC2 
show it satisfies (|05l) . So let H = 2rn' - 261 fo 
that H — c/r. But r^^H G L'^ad (since n' G L^^^) so c = 0, and hence H = i.e. p.45p holds. This implies 
immediately that n' e BC^ and n € BC^ with a bound ||n||Bci.i < c{\Sl\ + \\(3\\bc^ + IkllscO- 

Next define X2 G Hj:^'^^ by ([Qe]) : it follows that l(07|) holds, and so X2 is G^ for r > 0. To estabhsh 
X2 € BCl 
equation: 



Next for n which was defined to solve jS^Hl); we must 
q- l3, then (jOSll and ([QI]) imply that (riJ)' = so 



1,0 



it is necessary to analyze the behaviour at the origin. To achieve this, integrate up the fourth 



X2W 



'1 



mu-. 



)Xi-]2 



mn 



-ip 



s ds. 



The quantity in square brackets is continuous by assumption, so that lim,.^o r~^X2('') exists and is finite; it 
follows from (I3.47P that lim^-^o X2('') exists and is finite also, so that X2 "= BGJ'" as required. Finally p.46p 
then imphes that xi G BC^'°. □ 

Lemma 13 (Exponentially weighted bounds) Assume u € H^^j^{M.'^) solves 

{-A + M + V{r))u = f 

where f G Hj:°f 



), for some S < \JM. Assume further that V is continuous and limr^+00 V{r) = 0. 
Then u also satisfies the hound ||u||^{2,5} < c(||M||i/2 + ll/Uffio.'*}) for some c = c{6) > 0. 



Proof V = ru solves 



- v" + Mv + Vv = rf 



(3.49) 



11 



For any S > R + 1 > R let b{r) be a smooth function with b{r) — Oiir<Rorr>S+l, and b{r) = 1 if 
S > r > R+l. Multiply by e^^''h{r) and integrate, estimate |/ 25e^^''bvv' dr\ < J {e'^^^'S'^bv'^ /{I - e) + (1 - 
e)e^^^bv''^)dr and integrate by parts all other terms involving vv' . This leads to 



(1-e) 







dr. 



For any e > let R be such that sup^>^ \V{r)\ < e, and let 5^ < (a/ _ 2e)(l - e), so that \M-j^ + V\ > e 
for r > R. Now let 5 ^ +oo, to deduce 



R+l 



dr < c(e)(|lu|l|2 + / re2'5'-|/z;|dr). (3.50) 



Using the fact that < c(||m|1^i + jR+i^^^'' [v"^ +v''^)dr), (|330l) implies \\u\\h{i.s} < c{\\u\\h^ + 

||/||^{o,i}). To improve this to H^'^'^^ it is only necessary to multiply the equation (|3.49p by e^^b^ , square 
and integrate, to obtain 



R-\- 1 



since ||/||^{o,5} = 47r J Pe^^"^ pdr. But ||u||^{2,i} < c(||u||^,i + ^^j^^e'^^'^v"'^dr) since for an arbitrary 
second order derivative r^jViV^wp < cr^du''^ + r^^ju'p) so that in the region r > i? + 1 > 1 there holds 
r^jViVjUp < c(|u"p + r^lup + r^ju'p). This completes the proof. □ 
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